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Abstract 

It is well known that to determine a triangle up to congruence requires three measurements: three 
sides, two sides and the included angle, or one side and two angles. We consider various generalizations 
of this fact to two and three dimensions. In particular we consider the following question: given a convex 
polyhedron P, how many measurements are required to determine P up to congruence? 

We show that in general the answer is that the number of measurements required is equal to the 
number of edges of the polyhedron. However, for many polyhedra fewer measurements suffice; in the 
case of the unit cube we show that nine carefully chosen measurements are enough. 

We also prove a number of analogous results for planar polygons. In particular we describe a variety 
of quadrilaterals, including all rhombi and all rectangles, that can be determined up to congruence with 
only four measurements, and we prove the existence of n-gons requiring only n measurements. Finally, 
we show that one cannot do better: for any sequence of n distinct points in the plane one needs at least n 
measurements to determine it up to congruence. 

1 Introduction 

We discuss a class of problems about the congruence or similarity of three dimensional polyhedra. The basic 
format is the following: 

Problem 1.1. Given two polyhedra in R 3 which have the same combinatorial structure (e.g. both are 
hexahedra with four-sided faces), determine whether a given set of measurements is sufficient to ensure that 
the polyhedra are congruent or similar. 

We will make this more specific by specifying what sorts of measurements will be allowed. For example, 
in much of the paper, allowed measurements will include distances between pairs of vertices, angles between 
edges, angles between two intersecting face diagonals (possibly on different faces with a common vertex) or 
between a face diagonal and an edge, and dihedral angles (that is, angles between two adjoining faces). One 
motivation for these choices is given below. Sometimes we are more restrictive, for example, allowing only 
distance measurements. 

In two dimensions this was a fundamental question answered by Euclidean geometers, as (we hope) every 
student who takes geometry in high school learns. If the lengths of the corresponding sides of two triangles 
are equal, then the triangles are congruent. The SAS, ASA, and AAS theorems are equally well known. The 
extension to other shapes is not often discussed, but we will have some remarks about the planar problem as 
well. It is surprising to us that beyond the famous theorem of Cauchy discussed below, we have been unable 
to find much discussion of the problems we consider in the literature, though we think it is almost certain 
that they have been studied in the past. We would be appreciative if any reader can point us to relevant 
results. 

Our approach will usually be to look at the problem locally. If the two polyhedra arc almost congruent, 
and agree in a certain set of measurements, are they congruent? At first glance this looks like a basic 
question in what is known as rigidity theory, but a little thought shows that it is different. In rigidity 
theory, attention is paid to relative positions of vertices, viewing these as connected by inextensiblc rods 
which are hinged at their ends and so can rotate relative to each other, subject to constraints imposed by 



the overall structure of rods. In our problem there is the additional constraint that in any movement of the 
vertices, the combinatorial structure of the polyhedron cannot change. In particular, any vertices that were 
coplanar before the movement must be coplanar after the movement. This feature seems to us to produce 
an interesting area of study. 

Our original motivation for considering this problem came from a very practical question encountered 
by one of us (SPH). If one attempts to make solid wooden models of interesting polyhedra, using standard 
woodworking equipment, it is natural to want to check how accurate these models are. 1 As a mathematician 
one may be attracted first to the Platonic solids, and of these, the simplest to make appears to be the cube. 
(The regular tetrahedron looks harder, because non-right angles seem harder to cut accurately. ) 

It is possible to purchase lengths of wood with a square cross section, called "turning squares" because 
they are mostly used in lathes. To make a cube, all one has to do is use a saw to cut off the correct length 
piece from a turning square. Of course, one has to do so in a plane perpendicular to the planes of sides 
of the turning square. It is obvious that there are several degrees of freedom, meaning several ways to go 
wrong. The piece cut off could be the wrong length, or you could cut at the wrong angle, or perhaps the 
cross section wasn't square to begin with. So, you start measuring to see how well you have done. 

In this measurement, though, it seems reasonable to make some assumptions. The basic one of interest 
here is that the saw cuts off a planar slice. You also assume that this was true at the sawmill where the 
turning square was made. So you assume that you have a hexahedron - a polyhedron with six faces, all of 
which are quadrilaterals. Do you have a cube? At this point you are not asking a question addressed by 
standard rigidity theory. 

One's first impulse may be to measure all of the edges of the hexahedron, with the thought that if these 
are equal, then it is indeed a cube. This is quickly seen to be false, because the faces could be rhombi. 
Another intriguing possibility that we considered early on is that measuring the twelve face diagonals might 
suffice. However, we found some examples showing that this was not the case, and David Allwright [2] gave 
a short and elegant quaternion-based classification of all hexahedra with equal face diagonals. See also an 
earlier discussion of this problem in [13]. 2 Clearly some other configuration of measurements, perhaps 
including angles, is necessary. It did not take long to come up with several sets of 12 measurements which 
did the job, but a proof that this number was necessary eluded us. 

In our experience most people, even most mathematicians, who are presented with this problem do not see 
an answer immediately. Apparently the cube is harder than it looks, and so one would like a simpler problem 
to get some clues. The (regular) tetrahedron comes to mind, and so one asks how many measurements are 
required to establish that a polyhedron with four triangular faces is a tetrahedron. 

Now we turn to Cauchy's theorem. 

Theorem 1.2 (Cauchy, 1839). Two convex polyhedra with corresponding congruent and similarly situated 
faces have equal corresponding dihedral angles. 

If we measure the six edges of our triangular faced object, and find them equal, then we have established 
congruence of the faces of our object to those of a tetrahedron. Cauchy's theorem tells us that the dihedral 
angles are the same and this implies the desired congruence. 

For a tetrahedron this result is fairly obvious, but for the two other Platonic solids with triangular faces, 
namely the octahedron and the icosahedron, it is less so. Hence Cauchy's theorem is of practical value to 
the (extremely finicky) woodworker, and shows that for these objects, the number of edges is at least an 
upper bound on the number of measurements necessary to prove congruence. From now on we will denote 
the number of edges of our polyhedron by E, the number of vertices by V and the number of faces by F. 
We will only consider simply connected polyhedra, so that Euler's formula, V + F = E + 2, holds. 

It is not hard to give an example showing the necessity of convexity in Cauchy's result, but it is one where 
the two polyhedra being compared are, in some sense, far apart. It was not easy to determine if convexity 

lr rhe usual method of constructing a polyhedron is by folding a paper shell. 

2 The face diagonals of a hexahedron are equal when these diagonals form two congruent regular tetrahedra whose edges 
intersect in pairs. As it turns out, this arrangement is not unique. We have included Allwright's analysis of the possibilities as 
Appendix A in the posting of this paper to www.arXiv.org. 
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was necessary for local congruence. Can a nonconvex polyhedron with all faces triangular be perturbed 
smoothly through a family of noncongruent polyhedra while keeping the lengths of all edges constant? The 
answer is yes, as was proved in a famous and important paper by R. Connelly in 1979. [ ] 

Cauchy's result also gives us an upper bound for the number of measurements necessary to determine 
a unit cube: triangulate the cube by dividing each square face into a pair of triangles. Then we have a 
triangular-faced object with eighteen edges, and by Cauchy's theorem those eighteen edge measurements 
suffice to determine the cube up to congruence. 

However, we can do better. Start by considering a square. We can approach this algebraically by assuming 
that one vertex of the square is at (0, 0) in the plane. Without loss of generality we can also take one edge 
along the x-axis, going from (0,0) to (£i,0) for some x\ > 0. The remaining vertices are then (£2, £3) and 
(x^x*,) and this leads to the conclusion that to determine five unknowns, we need five equations, and so 
five measurements. For example, we could measure the four sides of the square and one vertex angle, or we 
could measure a diagonal instead of the angle. 

We then use this to study the cube. Five measurements show that one face is a square of a specific size. 
Only four more are needed to specify an adjacent face, because of the common edge, and the three more for 
one of the faces adjoining the first two. The requirement that this is a hexahedron then implies that we have 
determined the cube completely, with twelve measurements. This is a satisfying result because it shows that 
E measurements suffice for a cube as well as for the triangular faced Platonic solids. 

However, as remarked earlier, at this stage we have not proved the necessity of twelve measurements, 
only the sufficiency. One of the most surprising developments for us in this work was that in fact, twelve 
are not necessary. It is possible to determine a cube (including its size) with nine measurements of distances 
and face angles. The reason, even more surprisingly, is that only four measurements are needed to determine 
the congruence of a quadrilateral to a specific square, rather than five as seemed so obvious in the argument 
above. 

We will give the algorithm that determines a square in four measurements in the final section of the paper, 
which contains a number of remarks about congruence of polygons. For now, we proceed with developing 
a general method for polyhedra. This method will also handle similarity problems, where the shape of the 
polyhedron is specified up to a scale change. In determining similarity, only angle measurements are involved. 
As the reader might expect, in general E — 1 measurements suffice, with one additional length required to 
get congruence. 

2 E measurements suffice 

In this section we prove that for a convex polyhedron P with E edges, there is a set of E measurements 
that, at least locally, suffices to determine P up to congruence. 

We restrict to convex polyhedra mostly for reasons of convenience: many of the results below should be 
true in greater generality. (One problem with moving beyond convex polyhedra is determining exactly what 
the term 'polyhedron' should mean: for a recent attempt to give a general definition of the term 'nonconvex 
polyhedron', sec the beautiful paper [ ].) To avoid any ambiguity we begin with a precise definition of convex 
polyhedron. 

Definition 2.1. A closed half-space is a subset of R 3 of the form { (x, y, z) S R 3 | ax + by + cz + d > } 
with (a, b, c) ^ (0,0,0). A convex polyhedron is a subset P of R 3 which is bounded, does not lie in any plane, 
and can be expressed as an intersection of finitely many closed half-spaces. 

The vertices, edges and faces of a convex polyhedron P can be defined in terms of intersections of P with 
suitable closed half-spaces. For example, a face of P is a subset of P of the form P n H for some closed 
half-space H, such that P n H lies entirely within some plane but is not contained in any line. Edges and 
vertices can be defined similarly. 

The original problem refers to two polyhedra with the same 'combinatorial structure', so we give a notion 
of abstract polyhedron which isolates the combinatorial information embodied in a convex polyhedron. 
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Definition 2.2. The underlying abstract polyhedron of a convex polyhedron P is the triple (Vp, !Fp,Tp), 
where Vp is the set of vertices of P, Tp is the set of faces of P, and Ip C Vp x Tp is the incidence relation 
between vertices and faces; that is, (v, f) is in Ip if and only if the vertex v lies on the face /. 

Thus to say that two polyhedra P and Q have the same combinatorial structure is to say that their 
underlying abstract polyhedra are isomorphic; that is, there are bijections fly ■ Vp — > Vq and j3 F : Tp — > Tq 
that respect the incidence relation: (v, f) is in Tp if and only if (f3y(v), /3p(/)) is in Iq. Note that there 
is no need to record information about the edges; we leave it to the reader to verify that the edge data 
and incidence relations involving the edges can be recovered from the incidence structure (Vp,Tp,Tp). The 
cardinality of the set Tp is twice the number of edges of P, since 



and the latter sum counts each edge of P exactly twice. 

For the remainder of this section, we fix a convex polyhedron P and write V, E and F for the number 
of vertices, edges and faces of P, respectively. Let II = (V,T,1) be the underlying abstract polyhedron. We 
arc interested in determining which sets of measurements are sufficient to determine P up to congruence. 
A natural place to start is with a naive dimension count: how many degrees of freedom does one have in 
specifying a polyhedron with the same combinatorial structure as PI 

Definition 2.3. A realization of II = (V, J 7 , 1) is a pair of functions (av,ffy) where ay : V — > R 3 gives a 
point for each v in V, otjr; T — > {planes in R 3 } gives a plane for each / in T, and the point ay(v) lies on 
the plane ap(f) whenever (v, f) is in X. 

Given any convex polyhedron Q together with an isomorphism [3: (Vq,Tq,1q) = (V,T,X) of incidence 
structures we obtain a realization of II, by mapping each vertex of P to the position of the corresponding 
(under /?) vertex of Q and mapping each face of P to the plane containing the corresponding face of Q. In 
particular, P itself gives a realization of II, and when convenient we'll also use the letter P for this realization. 
Conversely, while not every realization of II comes from a convex polyhedron in this way, any realization 
of II that's sufficiently close to P in the natural topology for the space of realizations gives — for example by 
taking the convex hull of the image of ay — a convex polyhedron whose underlying abstract polyhedron can 
be identified with II. So the number of degrees of freedom is the dimension of the space of realizations of II 
in a neighborhood of P. 

Now we can count degrees of freedom. There are 3V degrees of freedom in specifying ay and 3P in 
specifying ujr. So if the |X| = IE 'vertex-on-face' conditions are independent in a suitable sense then the 
space of all realizations of II should have dimension 3V+3P— 2E or using Euler's formula — dimension E+6. 
We must also take the congruence group into account: we have three degrees of freedom available for 
translations, and a further three for rotations. Thus if we form the quotient of the space of realizations by 
the action of the congruence group, we expect this quotient to have dimension E. This suggests that E 
measurements should suffice to pin down P up to congruence. 

In the remainder of this section we show how to make the above naive dimension count rigorous, and 
how to identify specific sets of E measurements that suffice to determine congruence. The main ideas 
are: first, to use a combinatorial lemma (Lemma 2.7) to show that the linearizations of the vertex-on-face 
conditions are linearly independent at P, allowing us to use the inverse function theorem to show that the 
space of realizations really does have dimension E + 6 near P and to give an infinitesimal criterion for a 
set of measurements to be sufficient (Theorem 2.10), and second, to use an infinitesimal version of Cauchy's 
rigidity theorem to identify sufficient sets of measurements. 

The various measurements that we're interested in can be thought of as real-valued functions on the 
space of realizations of II (defined at least on a neighborhood of P) that are invariant under congruence. We 
single out one particular type of measurement: given two vertices v and w of P that lie on a common face, 
the face distance associated to v and w is the function that maps a realization Q = (ay,ajr) of II to the 
distance from a\>(v) to ay(w). In other words, it corresponds to the measurement of the distance between 
the vertices of Q corresponding to v and w. The main result of this section is the following theorem. 




(number of edges on /) 
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Theorem 2.4. Let P be a convex polyhedron with underlying abstract polyhedron {V,T,T). Then there is a 
set S of face distances of P such that (i) S has cardinality E, and (ii) locally near P, the set S completely 
determines P up to congruence in the following sense: there is a positive real number e such that for any 
convex polyhedron Q and isomorphism f3: (V,!F,T) = (Vq,J-q,Xq) of underlying abstract polyhedra, if 

1. each vertex v of P is within distance e of the corresponding vertex /3y(f) of Q, and 

2. m(Q) = m(P) for each measurement m in S , 
then Q is congruent to P. 

Rephrasing: 

Corollary 2.5. Let P be a convex polyhedron with E edges. Then there is a set of E measurements that is 
sufficient to determine P up to congruence amongst all nearby convex polyhedra with the same combinatorial 
structure as P. 

We'll prove this theorem as a corollary of Theorem 2.10 below, which gives conditions for a set of 
measurements to be sufficient. We first fix some notation. Choose numbcrings vi, . . . ,Wy and /i, . . . , fp 
of the vertices and faces of II, and write (xi(P),yi(P), Zi(P)) for the coordinates of vertex Vi of P. We 
translate P if necessary to ensure that no plane that contains a face of P passes through the origin. This 
allows us to give an equation for the plane containing fj in the form aj(P)x + bj(P)y + Cj(P)z — 1 for 
some nonzero triple of real numbers (aj(P),bj(P),Cj(P)); similarly, for any realization Q of II that's close 
enough to P the ith vertex of Q is a triple (xi(Q), yi(Q), Zi(Q) and the jth plane of Q can be described by an 
equation aj(Q)x+bj(Q)y+Cj(Q)z = 1. Hence the coordinate functions (x%, y±, z\,X2,y2, %2, ■ ■ ■ , cti, b%, c\, . . . ) 
give an embedding into R 3V+3F of some neighborhood of P in the space of realizations of II. 

For every pair fj) in T a realization Q should satisfy the 'vertex-on-face' condition 

a>j(Q)xi(Q) + b 3 (Q)y z (Q) + c,(Q)z 4 (Q) = 1. 

Let 4>ij be the function from R 3l/ + 3F to R defined by 

<j>i,j(xi,yi, zi,... ,oi,6i,ci, . . . ) = ajXi + bjyi + CjZ, L - 1, 

and let <j>\ ~R? V + 3F — > R 2B be the vector- valued function whose components are the 4>i,j as ( v ii fj) runs over 
all elements of T (in some fixed order). Then a vector in R' n + * u gives a realization of II if and only if it 
maps to the zero vector under <p. 

We next present a combinatorial lemma, Lemma 2.7, that appears as an essential component of many 
proofs of Steinitz's theorem, characterizing edge graphs of polyhedra. (See Lemma 2.3 of [ ], for example.) 
We give what we believe to be a new proof of this lemma. First, an observation that is an easy consequence 
of Euler's theorem. 

Lemma 2.6. Suppose that T is a planar bipartite graph of order r. Then there is an ordering rii, ri2, ■ ■ ■ , n r 
of the nodes of T such that each node rii is adjacent to at most three preceding nodes. 

Proof. We give a proof by induction on r. If r < 3 then any ordering will do. If r > 3 then we can apply a 
standard consequence of Euler's formula (see, for example, Theorem 16 of [3]), which states that the number 
of edges in a bipartite planar graph of order r > 3 is at most 2r — 4. If every node of T had degree at least 4 
then the total number of edges would be at least 2r, contradicting this result. Hence every nonempty planar 
bipartite graph has a node of degree at most 3; call this node n r . Now remove this node (and all incident 
edges), leaving again a bipartite planar graph. By the induction hypothesis, there is an ordering ni, . . . , rt r _i 
satisfying the conditions of the theorem, and then m, . . . , n r gives the required ordering. □ 

Lemma 2.7. Let P be a convex polyhedron. Consider the set V U T consisting of all vertices and all faces 
of P. Lt is possible to order the elements of this set such that every vertex or face in this set is incident with 
at most three earlier elements ofVUJ-. 



5 



Proof. We construct a graph T of order V + F as follows. T has one node for each vertex of G and one node 
for each face of G. Whenever a vertex v of P lies on a face / of P we introduce an edge of T connecting 
the nodes corresponding to v and /. Since P is convex, the graph T is planar; indeed, by choosing a point 
on each face of P, one can draw the graph T directly on the surface of P and then project onto the plane. 
(The graph T is known as the Levi graph of the incidence structure II = (V, J- ,T).) Now apply the preceding 
lemma to this graph. □ 

We now show that the functions fa j are independent in a neighborhood of P. Write Dcf>(P) for the 
derivative of <f> at P; as usual, we regard Dtfi(P) as a 2E-by-(3V + 3P) matrix with real entries. 

Lemma 2.8. The derivative D<fi(P) has rank 2E. 

In more abstract terms, this lemma implies that the space of all realizations of II is, in a neighborhood 
of P, a smooth manifold of dimension 3V + 3P — 2E = E + 6. 

Proof. We prove that there are no nontrivial linear relations on the 2E rows of Dcf)(P). To illustrate 
the argument, suppose that the vertex v\ lies on the first three faces and no others. Placing the rows 
corresponding to fa\, 4>\ t 2 and ^1.3 first, and writing simply x\ for x\(P) and similarly for the other 
coordinates, the matrix D<fi(P) has the following structure. 
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Here the vertical double bar separates the derivatives for the vertex coordinates from those for the face 
coordinates. Since the faces /1, f2 and fs cannot contain a common line, the 3-by-3 submatrix in the top 
left corner is nonsingular. Since v\ lies on no other faces, all other entries in the first three columns are zero. 
Thus any nontrivial linear relation of the rows cannot involve the first three rows. So D(f>(P) has full rank 
(that is, rank equal to the number of rows) if and only if the matrix obtained from D(j){P) by deleting the 
first three rows has full rank — that is, rank 2E — 3. Extrapolating from the above, given any vertex that lies 
on exactly three faces, the three rows corresponding to that vertex may be removed from the matrix D(f>(P), 
and the new matrix has full rank if and only if D<j)(P) does. The dual statement is also true: exchanging 
the roles of vertex and face and using the fact that no three vertices of P are collinear we see that for any 
triangular face / we may remove the three rows corresponding to / from D<j)(P), and again the resulting 
matrix has full rank if and only if the D<j)(P) does. Applying this idea inductively, if every vertex of P lies on 
exactly three faces (as in for example the regular tetrahedron, cube or dodecahedron), or dually if every face 
of P is triangular (as in for example the tetrahedron, octahedron or icosahedron) then the lemma is proved. 
For the general case, we choose an ordering of the faces and vertices as in Lemma 2.7. Then, starting from 
the top end of this ordering, we remove faces and vertices from the list one- by-one, removing corresponding 
rows of D(f)(P) at the same time. At each removal, the new matrix has full rank if and only if the old one 
does. But after removing all faces and vertices we're left with a 0-by-2P matrix, which certainly has rank 0. 
So DfaP) has rank 2E. □ 

We now prove a general criterion for a set of measurements to be sufficient. Given the previous lemma, 
this criterion is essentially a direct consequence of the inverse function theorem. 

Definition 2.9. A measurement for P is a smooth function m defined on an open neighborhood of P in 
the space of realizations of II, such that m is invariant under rotations and translations. 

Given any such measurement m, it follows from Lemma 2.8 that we can extend m to a smooth function on 
a neighborhood of P in R 3V+3F . Then the derivative Dm(P) is a row vector of length 3V + 3P, well-defined 
up to a linear combination of the rows in D(j){P). 
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Theorem 2.10. Let S be a finite set of measurements for H near P. Let ip: R 3V+SF _^ pJS| ^ e ^ e vec f or . 
valued function obtained by combining the measurements in S, and write D%p{P) for its derivative at P, an 
\S\-by-(3V + 3F) matrix whose rows are the derivatives Dm(P) for m in S. Then the matrix 

has rank at most 3E, and if it has rank exactly 3E then the measurements in S are sufficient to determine 
congruence: that is, for any realization Q ofH, sufficiently close to P, if m(Q) = m(P) for all m in S then 
Q is congruent to P. 

Proof. Let Q(t) be any smooth one-dimensional family of realizations of II such that Q(0) = P and Q(t) 
is congruent to P for all t. Since each Q(t) is a valid realization, (j>(Q(t)) = for all t. Differentiating and 
applying the chain rule at t = gives the matrix equation D<p(P)Q' (0) = where Q'(0) is thought of as a 
column vector of length 3V + 3F. The same argument applies to the map ip: since Q(t) is congruent to P for 
all t, ip(Q(t)) = ip(P) is constant and Dip(P)Q' (0) = 0. We apply this argument first to the three families 
where Q(t) is P translated t units along the rc-axis, y-axis, or z-axis respectively, and second when Q(t) 
is P rotated by t radians around the x-axis, the y-axis and the z-axis. This gives 6 column vectors that 
are annihilated by both D(p(P) and Dip(P). Writing G for the (3V + 3F)-by-6 matrix obtained from these 
column vectors, we have the matrix equation 

It's straightforward to compute G directly; we leave it to the reader to check that the transpose of G is 
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For the final part of this argument, we introduce a notion of normalization on the space of realizations of II. 
We'll say that a realization Q of II is normalized if v\{Q) = V\(P), the vector from V\(Q) to v 2 (Q) is in the 
direction of the positive x-axis, and Vi(Q), v 2 {Q) and v 3 (Q) all lie in a plane parallel to the xy-plane, with 
f 3 (Q) lying in the positive y-direction from V\(P) and ^(P). In terms of coordinates we require that X\{P) = 
Xi(Q) < x 2 (Q), yi(P) = yi(Q) = y 2 {Q) < ya{Q) and Zl (P) = z^Q) = z 2 (Q) = z 3 {Q). Clearly every 
realization Q of II is congruent to a unique normalized realization, which we'll refer to as the normalization 
of Q. Note that the normalization operation is a continuous map on a neighborhood of P. Without loss of 
generality, rotating P around the origin if necessary, we may assume that P itself is normalized. The condition 
that a realization Q be normalized gives six more conditions on the coordinates of Q, corresponding to six 
extra functions Xi ; • • ■ ,X6, which we use to augment the function (4>,tp) : R 3V+3i? — > R 2£ +I s l to a function 
0,V>,x): R 3y+3F -» R 2£ +ls|+6. xhcse six f unc tions are simply X i(Q) = MQ) ~ M p ), X2(Q) = Vi(Q) - 
Vi(P), X3(Q) = z x {Q) - zt(P), X i(Q) - V2(Q) - 2/1 OP), X6(Q) = Z2(Q) - zi(P) and Xe (Q) = z 3 (Q) - Zl (P). 
Claim 2.11. The 6-by-6 matrix Dx(P)G is invcrtible. 

Proof. The matrix for G was given earlier; the product D\{P)G is easily verified to be 
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which has nonzero determinant (2/3 — y\){x2 — x\) 2 . 



□ 



As a corollary, the columns of G are linearly independent, which proves that the matrix in the statement of 
the theorem has rank at most 3E. Similarly, the rows of D\(P) must be linearly independent, and moreover 
no nontrivial linear combination of those rows is a linear combination of the rows of Dip(P). Hence if DiJj(P) 
has rank exactly 3E then the augmented matrix 



has rank 3-E+6. Hence the map (0, ip, x) has injective derivative at P, and so by the inverse function theorem 
the map (0, tp, x) itself is injective on a neighborhood of P in R 3V+3F . Now suppose that Q is a polyhedron 
as in the statement of the theorem. Let R be the normalization of Q. Then (j>(R) — <j>(Q) — 4>{P) = 0, 
'tp(R) — ip(Q) — ip(P), and x(R) = x(P)- So if Q is sufficiently close to P, then by continuity of the 
normalization map R is close to P and hence R — P by the inverse function theorem. So Q is congruent 
to R = P as required. □ 

Definition 2.12. Call a set S of measurements sufficient for P if the conditions of the above theorem apply: 
that is, the matrix D((p,ijf)(P) has rank 3E. 

Corollary 2.13. Given a sufficient set S of measurements, there's a subset of S of size E that's also 
sufficient. 

Proof. Since the matrix D(4>,i/j)(P) has rank 3E by assumption, and the 2E rows coming from D<fi(P) are 
all linearly independent by Lemma 2.8, we can find E rows corresponding to measurements in 5' such that 
Dcf>(P) together with those E rows has rank 3E. □ 

The final ingredient that we need for the proof of Theorem 2.4 is the infinitesimal version of Cauchy's 
rigidity theorem, originally due to Dehn, and later given a simpler proof by Alexandrov. We phrase it in 
terms of the notation and definitions above. 

Theorem 2.14. Let P be a convex polyhedron, and suppose that Q(t) is a continuous family of polyhedra 
specializing to P at t = 0. Suppose that the real-valued function mo Q is stationary (that is, its derivative 
vanishes) at t — for each face distance m. Then Q'(0) is in the span of the columns of G above. 

Proof. See Chapter 10, section 1 of [1]. See also section 5 of [7] for a short self-contained version of Alexan- 
drov's proof. □ 

Corollary 2.15. The set of all face distances is sufficient. 

Proof. Let S be the collection of all face distances. Then the matrix D(4>,ip)(P) has rank 3E. Now Theo- 
rem 2.14 implies that for any column vector v such that D(4>,ip)(P)v — 0, v is in the span of the columns 
of G. Hence the kernel of the map D(<f), ip)(P) has dimension exactly 6 and so by the rank- nullity theorem 
together with Euler's formula the rank of D(4>,ip)(P) is 3V + 3F — 6 = 3E. □ 

Now Theorem 2.4 follows from Corollary 2.15 together with Corollary 2.13. 

Finding an explicit sufficient set of E face distances is now a simple matter of turning the proof of 2.13 
into a constructive algorithm. First compute the matrices Dtfi(P) and Dip(P) (the latter corresponding to 
the set S of all face distances). Initialize a variable T to the empty set. Now iterate through the rows of 
Dip(P): for each row, if that row is a linear combination of the rows of D(j>(P) and the rows of Dip(P) 
corresponding to measurements already in T, discard it. Otherwise, add the corresponding measurement 
to T. Eventually, T will be a sufficient set of E face distances. 

A computer program has been written in the language Python to implement this algorithm. This program 
is attached as Appendix B of the online version of this paper, available at www.arXiv.org. We hope that 



( 



Dct>{P) 
Di/}(P) 
D X (P) 



) 
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the comments within the program are sufficient explanation for those who wish to try it, and we thank Eric 
Korman for a number of these comments. 

The algorithm works more generally. Given any sufficient set of measurements, which may include angles, 
it will extract a subset of E measurements which is sufficient. It will also find a set of angle measurements 
which determines a polyhedron up to similarity. As an example we consider a similarity calculation for a 
dodecahedron. Our allowed measurements will be the set of angles formed by pairs of lines from a vertex 
to two other points on a face containing that vertex. In Figure 2 we show a set of 29 such angles which the 
program determined to characterize a dodecahedron up to similarity. 



There is no restriction of the method to Platonic solids. Data for a number of other examples can be 
found in the program listing. Among these examples are several which are non-convex. The program appears 
to give a result for many of these, but we have not extended the theory beyond the convex case. 

2.1 Related work 

The results and ideas of section 2 are, for the most part, not new. The idea of an abstract polyhedron 
represented by an incidence structure, and its realizations in R 3 , appears in section 2 of [14]. In Corollary 15 
of that paper, Whiteley proves that the space of realizations of a 'spherical' incidence structure (equivalent 
to an incidence structure arising from a convex polyhedron) has dimension E. The essential combinatorial 
content of the proof of Lemma 2.8 is often referred to as 'Steinitz's lemma', and a variety of proofs appear 
in the literature ([10], [8]); we believe that the proof above is new. 

3 How many measurements are enough? 

Theorem 2.4 provides an upper bound for the number of distance/ angle measurements needed to describe a 
polyhedron with given combinatorial structure. But it turns out that many interesting polyhedra can be de- 
scribed with fewer measurements. In particular, a cube can be determined by 9 distance/angle measurements 
instead of 12. Also, 10 distance measurements (no angles used) suffice. 

This phenomenon appears already in dimension two. Much of this section deals with polygons and sets of 
points in the plane. In this simpler setting one can often find precisely the smallest number of measurements 
needed. Extending these results to polyhedra, with or without fixing their combinatorial structure, certainly 
deserves further study. 




Figure 2. A set of 29 angles 
which locally determine a 
dodecahedron up to 
similarity. 
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So, we turn our attention to polygons. Unlike polyhedra, their combinatorial structure is relatively simple, 
especially if we restrict to the case of convex polygons. One can argue about the kinds of measurements that 
should be allowed, but the following three seem the most natural to us. 

Definition 3.1. Suppose that A\A 2 ■ ■ ■ A n is a convex polygon. (More generally, suppose {Ai, A2, A n } 
is a sequence of distinct points on the plane). Then the following are called simple measurements: 

1. distances l-AjA^I , i ^ j 

2. angles ZAiAjAk for i,j, k distinct 

3. 'diagonal angles' between AiAj and A^Ai. 

Other quantities might also be considered, like the distance from Ai to AjA^, but in practice these would 
require several measurements. 

It is natural to ask how many simple measurements arc needed to determine a polygon up to isometry. 
In the case of triangle the answer is 3. There are a couple of different ways to do it, which are very important 
for the classical Euclidean Geometry (SSS, SAS, SSA, ASA, AAS). It is appropriate to note that two sides 
of a triangle and an angle adjacent to one of the sides do not always determine the triangle uniquely. But 
they do determine it locally (as in part ii of Theorem 2.4). 

It is easy to see that for any n— gon P = A\ ■ ■ ■ A n the following (2n — 3) measurements suffice: the 
Ti—l distances |AiA 2 |, |A 2 A 3 |, . . . , |A„_x^4 n |, together with the n — 2 angles ZAiA 2 A 3 , . . . , ZA n _ 2 A n _iA n . 
Observe also that instead of using distances and angles, one can use only distances, by substituting 
for ZA i - 1 A i A i+1 . 

The following theorem implies that for most polygons one can not get away with fewer measurements. 
To make this rigorous, we will assume that A\ = (0,0) , and that A 2 — (x 2 ,0) for some x 2 > 0. It then 
becomes obvious that to each polygon we can associate a point in R 2 ' 1-3 , corresponding to the undetermined 
coordinates x 2 , x 3 ,y 3 ,..., x n , y n with (xi,yi) ^ (xj,Vj) for i ± j. 

Theorem 3.2. Denote by S the set of all points in R 2,l ~ 3 obtained by the procedure above from those 
polygons that can be determined up to isometry by fewer than 2n — 3 measurements. Then S has measure 
zero. 

Proof. Observe that any set of 2n — 3 specific measurements is a smooth map from R 2 "~ 3 to itself. There is 
only a finite number of such maps, with measurements chosen from the the types 1,2 and 3. At a non-critical 
point, the map has an inverse. The result is then a consequence of Sard's theorem (Chapter 2, Theorem 8 
of [12]). □ 

The above theorem is clearly not new, and we claim no originality for its statement or proof. For n = 3 it 
implies that for a generic triangle one needs at least three simple measurements. This is true for any triangle, 
because two measurements are clearly not enough. For 71 = 4 the theorem implies that a generic convex 
polygon requires 5 measurements. One can be tempted to believe that no quadrilateral can be described by 
fewer than 5 measurements. In fact, the first reaction of most mathematicians, seems to be that if a general 
case requires 5 measurements, all special cases do so as well. The following example seems to confirm this 
intuition. 

Example 3.3. Suppose P — AiA 2 A 3 A4 is a quadrilateral. Suppose ZA 2 A\A^ — a%, ZA3A4A1 = a 2 , 
\AiA 2 \ = di, \A 3 Ai\ = d 2 , and ZA\A 2 A 3 = a 3 . For most polygons P these five measurements are sufficient 
to determine P up to congruence. However for some P they are not sufficient. For example, if P is a square, 
there arc infinitely many polygons with the same measurements: all rectangles with the same |j4i |- 

The above example suggests that the special polygons require at least as many measurements as the 
generic ones. So (2n — 3) should be the smallest number of simple measurements required for any n— gon. 
However, despite the above example, this is very far from the truth: many interesting n— gons can be 
determined by fewer than (2n — 3) simple measurements. 
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The simplest example in this direction is not usually considered a polygon. Suppose A\A 2 A 3 A 4 is a set 
of 4 points on the plane that lie on the same line, in the natural order. Then the 4 distances |.AiAi|, 1^4x^2 1 , 
1^42^3 1^3^41 determine the configuration up to the isometry. This means that any subset of four points 
on the plane with the same corresponding measurements is congruent to A 1 A 2 A 3 A i . The proof of this is, of 
course, a direct application of the triangle inequality: 

\AtA4\ < \AiAsl + \A 3 A 4 \ < \AxA 2 \ + \AxA 2 \ + \A 2 A 3 \ + \A 3 A 4 \ 

One can think of A X A 2 , A 2 A 3 , A 3 A±, and A X A 4 as of rods of fixed length. When \A\A4\ = \AiA 2 \ + \A 1 A 2 \ + 
\A 2 A 3 \ + I v4_3 v4_4 1 , the configuration allows no room for wiggling, all rods must be lined up. This degenerate 
example easily generalizes to n points and n distance measurements. More surprisingly, it generalizes to 
some convex polygons. 

Definition 3.4. A polygon is called exceptional if it can be described locally by fewer than 2n — 3 measure- 
ments. 

Obviously, no triangles are exceptional, in the above sense. But already for the quadrilaterals, there exist 
some exceptional ones, that could be defined by four rather than five measurements. The biggest surprise 
for us was the following observation. 

Proposition 3.5. All squares are exceptional! 

Specifically, for a square ABCD the following four measurements distinguish it among all quadrilaterals: 

\AB\,\AC\,\AD\,ZBCD. 

Proof. Suppose A'B'C'D' is another quadrilateral with \A'B'\ = \AB\, \A'C'\ = \AC\, \A'D'\ = \AD\, and 
ZB'C'D' = ABCD. If \AB\ = d, this means that \A'B'\ = \A'D'\ = d, \A'C'\ = dV2, and ZB'C'D' = f . We 
claim that the quadrilateral A'B'C'D' is congruent to the square ABCD. Indeed, consider A' fixed. Then 
B' and D' lie on the circle of radius d centered at A' , while C lies on the circle of radius d^/2 centered at A' . 
Then ^ is the maximal possible value of the angle B'C'D' and it is achieved when ZA'B'C = ZA'D'C = |, 
making A'B'C'D' a square, with the side \A'B'\ = d. □ 

The idea of showing that a particular polygon occurs when some angle or length is maximized within 
given constraints, and that there is only one such maximum, is at the heart of all of the examples in this 
section. It is related to the notion of second order rigidity of tensegrity networks, see [6] and [11]. 

One immediate generalization of this construction is the following 3-parameter family of exceptional 
quadrilaterals. 

Proposition 3.6. Suppose that ABCD is a convex quadrilateral with ZABC = ZCDA — |. Then ABCD 
is exceptional. It is determined up to congruence by the following four measurements: \AB\ , \AD\ , \AC\ , 
and ZBCD. 

Proof. Consider an arbitrary quadrilateral ABCD with the same four measurements as the quadrilateral 
that we are aiming for (note that we are not assuming that the ZABC = ZCDA = |). We can consider 
the points A and C fixed on the plane. Then the points B and D lie on two fixed circles around A, with 
radius \AB\ and radius \AD\. Among all such pairs of points on this circle, on opposite sides of AC, the 
maximum possible value of ZBCD is obtained for only one choice of B and D. This is the choice which 
makes CB and CD tangent to the circle at B and D, and so ZABC and ZADC are right angles, and 
the quadrilateral ABCD is congruent to the one we are aiming for. This family of quadrilaterals includes 
all rectangles. In a sense it is the biggest possible: one cannot hope for a four-parameter family requiring 
just four measurements. Another such family is given below. Note that it includes all rhombi that are not 
squares. 

Proposition 3.7. Suppose ABCD is a convex quadrilateral, and for given B and D, and given acute 
angles 81, 9 2 , choose A and C so that ZDAB = B\, ZDCB = 9 2 , and \AC\ is as large as possible. Then this 
determines a unique quadrilateral, up to congruence, and this quadrilateral has the property that \AB\ = \AD\ 
and \CB\ = \CD\ . 
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□ 



We leave the proof to the reader. It implies that for the set of quadrilaterals ABCD such that AC is a 
perpendicular bisector of BD, and the angles ZDAB and ZDCB are acute, the measurements \BD\ , \ AC\ , 
ZDAB and ZDCB are sufficient to determine ABCD. 

As David Allwright pointed out to the authors, one can further extend this example to the situation 
when ZDAB + ZDCB < f . 

The existence of so many exceptional quadrilaterals suggested that for bigger n it may be possible to go 
well below the (2n — 3) measurements. This is indeed correct, for every n there are polygons that can be 
defined by just n measurements. 

Proposition 3.8. Suppose that A\A 2 ■ ■ ■ A n is a convex polygon, with 

ZA X A 2 A Z = ZA X A Z A± = ■■■ = ZA x A k A k+1 = ■■■ = ZA x A n - X A n = | 

(Note that many such polygons exist for every n). Then A\A 2 ■ • • A n is exceptional, moreover the distances 
\AiA 2 \, \A\A n \ and the angles ZA k A k+ \A\, 2 < k < n — 1 determine the polygon. 

Proof. Suppose ZA k A k +\A\ = a k , 2 < k < (n — 1). Then because of the Law of Sines for the triangles 
AAiA k Ak+i, we obtain the following sequence of inequalities: 

\A lAn \<^ A ^< |AA - 2 ' <■■■< ^ 

sina„_i sina„_i • sina„_ 2 sina„_i sina 2 

Equality is achieved if and only if all angles AiA k A k +i are right angles, which implies the result. 
One can ask whether an even smaller number of measurements might work for some very special polygons. 
The following theorem shows that in a very strong sense the answer is negative. 

Theorem 3.9. For any sequence of distinct points A\, A2, . . . , A n on the plane, one needs at least n distance 
/ angle / diagonal angle measurements to determine it up to plane isometry. 

Proof. At least one distance measurement is needed, and we can assume that it is |j4i A2I. We assume that 
AiA 2 is fixed, so the positions of the other n — 2 points determine the set up to isometry. We identify the 
ordered set of coordinates of these points with a point in R 2 (" -2 ) = R 2 ™~ 4 . The set of all sequences of 
distinct points then corresponds to an open set U C R 2 ™~ 4 . Each measurement of type 1, 2, or 3 determines 
a smooth submanifold V C U. The following observation is the main idea of the proof. 

Lemma 3.10. Suppose that x € V. Then there exists an affine subspace W o/R 2n ~ 4 , of dimension at least 
2n — 6, which contains x and is such that for some open ball B containing x we have W C\ B C V . 

□ 

Proof. The proof of this lemma involves several different cases, depending on the kind of measurements used 
and whether A\ and/or A 2 arc involved. We give some examples, the other cases being similar. First suppose 
that the polygon is the unit square, with vertices at A\ — (0, 0) , A 2 = (1,0), A 3 = (1,1), and A4 — (0, 1) 
Suppose that the measurement is the angle ZA 2 AiA 3 . With Ai and A 2 fixed, we are free to move A 3 along 
the line y = x, and A4 arbitrarily. In this case, then, W could be three dimensional, one more than promised 
by the Lemma. Second, again with the unit square, suppose that the measurement is the distance from A3 
to A 4 . In this case, we can move A 3 and A 4 the same distance along parallel lines. Thus, 

W = {(l + c,l + d, c,l + d)} 

which is the desired two-dimensional affine space. Finally, suppose that n > 5 and that the measurement is 
the angle between A\A k and AiA m where 1,2, k, l,m are distinct. In this case we can move Ai arbitrarily, 
giving two free parameters, and we can move A m so that AiA m remains parallel to the original line containing 
these points. Then the length |j4jA m | can be changed. This gives a third parameter. Then the length A^A k 
can be changed, giving a fourth, and the remaining n— 5 points can be moved, giving 2n— 10 more dimensions. 
In this case, W is 2n — 6 dimensional. We leave other cases to the reader. □ 
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Continuing the proof of Theorem 3.9, we first show that k < n — 3 measurements is insufficient to 
determine the points A±,... ,A n . Recall that one measurement, l^i^l, was already used. Suppose that a 
configuration x € R 2 ™ -4 belongs to all the varieties V\,. . . ,Vk- We will prove that x is not an isolated point 
in V = n^Vi. 

Proof. Denote the affine subspaces obtained from Lemma 3.10 corresponding to V$ and x by Wj. Let 
W = W\ ("I • ■ • fl Wfc. Since each Wi has dimension at least 2n — 6, its codimension in R 2 ™~ 4 is less than or 
equal to 2. 3 Hence 

codimVF < codim W\ H h codim W fc < 2A; < 2 (n - 3) < 2n - 4, 

so IF must have dimension at least 2. Since W is contained in all of the varieties V\, .., Vfe, there is a two- 
dimensional affine subspace of points containing x and satisfying all of the measurements. A neighborhood of 
x in this subspace is contained in U, showing that x is not isolated. The case k = n — 2 is trickier, because the 
dimensional count does not work in such a simple way. In this case, it is possible that W — WilT-nWfe = {x}. 
However, W' — W% H ■ • • fl Wk must have dimension at least 2. We now consider the space W\ in more detail. 
Denote the measurement defining V\ by fx : R 2n ~ 4 — ► R. The linearization of fi at x is its derivative, 
Dfi (x) . For measurements of types 1, 2, or 3 above it is not hard to see that Dfi (x) ^ 0. If X is the 
null space of the 1 x (2n — 4) matrix Dfi (x) , then the tangent space of V\ at x can be defined as the affine 
subspace T x (Vi) — x + X. This has dimension 2n — 5. It is possible that W\ — T x (Vi) . This is the case 
in the first example in the proof of the Lemma above. Then, dim W\ — 2n — 5. Even if dim (W) = 2, we 
must have dim(14 7 ) > 1, and since W n U C V, x is not isolated in V. However in the second example in 
the Lemma above, (a distance measurement for the square), dim Wi — 2 < dimT^ (Vi) = 3. Suppose we add 
as the last two measurement the distance from A4 to Ai. Then one can easily work out the spaces exactly, 
and show that V\ fl W, when projected onto the first three coordinates (#3,2/3,2:4), is the intersection of a 
vertical cylinder with the x^,X4 plane. Again, x is not isolated in V. We have left to consider the general 
case, where dim W\ = 2n — 6. (We are assuming that in the lemma we always choose the maximal W.) Since 
dim T x (V\) = 2n — 5 and dim W' = 2, W' D T x (Vi) must be one dimensional. We wish to show that x is not 
isolated in V\ D W . We consider the map gi = fi\w'- If Dgi 0, then the implicit function theorem implies 
that is not an isolated zero of gi in W , proving the theorem. To show that Dgi ^ we can use the chain 
rule. Let i : W — > R 2 ™ -4 be the identity on W' . Then g\ = f\o i and so by the chain rule, Dgi = Dfi o i. 
If Dgi = 0, then (Dfi) \ w = 0, which implies that W C T x (Vi) . Since dimVF' = 2, this contradicts the 
earlier assertion that W' n T x (Vi) is one dimensional, completing the proof of the theorem. □ 

Similar ideas can be used to construct other interesting examples of exceptional polygons and polyhedra. 
The following are worth mentioning. 

1. There exist tetrahedra determined by just 5 measurements, instead of the generic 6. In particular, 
if measurement of dihedral angles is permitted, the regular tetrahedron can be determined using 5 
measurements. It seems unlikely that 4 would ever work. 

2. There exist 5-vertex convex polyhedra that are characterized by just 5 measurements. Moreover, four 
of the vertices are on the same plane, but, unlike in the beginning of the paper, we do not have to 
assume this a priori! To construct such an example, we start with an exceptional quadrilateral ABCD 
as in Proposition Proposition 3.6, with \AD\ < \AB\. We then add a vertex E outside of the plane 
ABCD, so that ZADE — ZAEB — J. There are obviously many such polyhedra. Now notice that 
the distances \AD\, \AC\, and angles ZAED, ZABE, ZBCD completely determine the configuration. 

3. As pointed out earlier, using four measurements for a square, one can determine the cube with just 9 
distance / angle measurements. Interestingly, one can also use 10 distance measurements for a cube. If 

3 If V is an affine subspace of R m , then it is of the form x + X, where X is a subspace of R. m . If the orthogonal complement 
of X in R, m is Y, then the codimension of V is the dimension of Y. 
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A\B\C\D\ is its base and A2B2C2D2 is a parallel face, then the six distances between A\, B\, D\ and 
A2 completely fix the comer. Then IA1C2I, I-B2C2I, I-D2C2I, and IC1C2I determine the cube, because 

|^iC 2 | 2 < \B 2 C 2 \ 2 + \D 2 C 2 \ 2 + |CiC 2 | 2 , 

with the equality only when the segments B2C2, D2C2, and C1C2 are perpendicular to the correspond- 
ing faces of the tetrahedron A 1 i? 1 Z? 1 A 2 . 

The last example shows that the problem is not totally trivial even when we only use the distance 
measurements. It is natural to ask for the smallest number of distance measurements needed to fix a non- 
degenerate n— gon in the plane. We have the following theorem. 

Theorem 3.11. Suppose A\, A2, ■ ■ . , A n are points in the plane, with no three of them on the same line. 
Then one needs at least min(2n — 3, 4p) distance measurements to determine it up to a plane isometry. 

Proof. Denote by k(n) the smallest number of distance measurements that allows us to fix some non- 
degenerate configuration as above. We need to show that k(n) > min(2n — 3, ^p). We use induction on n. 
So suppose that the statement is true for all sets of fewer than n points, and that there is a non-degenerate 
configuration of n points, Ai, A2, . . . ,A n that is fixed by k measurements where k < 2n — 3 and k < 4p. 
Because k < 4p, there is some point Aj, which is being used in less than 3 measurements. If it is only used 
in 1 measurement, the configuration is obviously not fixed. So we can assume that it is used in two distance 
measurements, |j4»j4i| and |AjA2|. Because A\, A2, and Ai are not on the same line, the circle with radius 
|Ajj4i| around Ai and the circle with radius | ^4^ ^4 2 1 around A 2 intersect transversally at Ai. Now remove Ai 
and these two measurements. By the induction assumption, the remaining system of (n — 1) points admits 
a small perturbation. This small perturbation leads to a small perturbation of the original system. Q. E. D. 

For n < 7 the above bound coincides with the upper bound (2n — 3) for the minimal number of measure- 
ments (see note after Theorem 1). One can show that for n > 8, is the optimal bound, with the regular 
octagon being the simplest 'distance-exceptional' polygon. In fact there are at least two different ways to 
define the regular octagon with 12 measurements. 

Example 3.12. Suppose AiA 2 A 3 . . .Ag is a regular octagon. Its 8 sides and the diagonals |Ai^4 5 |, | ^4 8 ^4 6 1 , 
I^^UI, and, finally, IA3A7I determine it among all octagons. A proof is elementary but too messy to include 
here, its main idea is that IA3A7I is the biggest possible with all other distances fixed. In fact, one can show 
that the distance between the midpoints of A1A5 and AqAs is maximal when ZA^AiAg = ZAqA^Ai = 

Example 3.13. Suppose A1A2A3 . . . A% is a regular octagon. Its 8 sides and the 4 long diagonals determine 
it among all octagons. This follows from [ ], Corollary 1 to Theorem 5. In fact, this example can be 
generalized to any regular (2n)— gon: a tensegrity structure with cables for the edges and struts for the long 
diagonals has a proper stress due to its symmetry and thus is rigid. This provides an optimal bound for the 
case of even number of points, and one can get the result for the odd number of points just by adding one 
point at a fixed distance to two vertices of the regular (2n)— gon. We should note that tensegrity networks 
and their generalizations have been extensively studied, see for example [6]. 

There are many open questions in this area, some of which could be relatively easy to answer. For example, 
we do not know whether either of the results for the cube (9 measurements, or 10 distance measurements) 
is best possible. Also, it is relatively easy to determine a regular hexagon by 7 measurements, but it is not 
known if 6 suffice. 

One can also try to extend the general results of this section to the three-dimensional case. Methods 
of Theorem 3.11 produce a lower estimate of (2n) for the number of distance measurements needed to 
determine the set of n points, no four of which lie on the same plane. One should also note that many of the 
exceptional polygons that we have constructed, for instance the square, are unique even when considered in 
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the three-dimensional space: the measurements guarantee their planarity. The question of determining the 
smallest number of measurements, using distances and angles, with the planarity conditions for the faces, 
seems to be both the hardest and the most interesting. But even without the planarity conditions or without 
the angles the answer is not known. 
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A Quadrilateral-faced hexahedrons with all face-diagonals equal 



by David Allwright 

In this appendix we classify hexahedrons with quadrilateral faces with all face-diagonals of length d. We 
shall show that such a hexahedron lies in cither (a) a 1-paramctcr family with dihedral symmetry of order 6 
and all the faces congruent, or (b) a 2-parameter family with a plane of symmetry and 2 congruent opposite 
faces joined by 4 symmetric trapezoids. The cube is a special case of both families and in fact has the 
maximum volume and the maximum surface area. We first construct the two families and then show that 
the classification is complete. 

To construct the first family let ABCD be a regular tetrahedron with base BCD in the xy-plane, vertex 
A on the z-axis, and edges of length d. Let I be a line parallel to the base and passing through the z-axis, and 
choose I such that the rotation of the line segment AB through ir about / intersects the line segment CD. 
This is a single constraint on a 2-parameter family of lines so there is a 1-parameter family of such lines 
Then let A', B', C, D' be the rotations of A, B, C, D through n about I. By choice of I, the points A'CB'D 
are coplanar and form a convex quadrilateral with both diagonals of length d. But rotation through 7r about 
I and rotation through 27r/3 about the z-axis together generate a dihedral group of order 6 permuting these 
8 points. Then the images of the quadrilateral A'CB'D under that group are the 6 congruent faces of a 
hexahedron with all face-diagonals of length d. 

To construct the second family, let ABCD be a regular tetrahedron with edge length d, and let P be 
a plane with AB on one side and CD on the other, and such that the reflection of AB in P intersects 
CD. This is a single constraint on a 3-paramctcr family of planes, so there is a 2-parametcr family of 
such planes P. Then let A', B' , C , D' be the reflections of A,B,C,D in P. By choice of P, the points 
A'CB'D form a quadrilateral with both diagonals of length d. Its reflection in P has the same property. 
The points AA'C'C then have AA' parallel to CC (both perpendicular to P) so they are coplanar and form 
a symmetric trapezoid with both diagonals d. The same holds for the other faces that cut P, so again we 
have a hexahedron with the required property. 

To complete the classification we now show than any hexahedron with all face-diagonals equal must lie 
in one of these families. So, suppose we have such a hexahedron, let ao be one vertex and let ai, a 2 , a 3 be 
the other ends of the diagonals of the 3 faces meeting at a . Then a 1; a 2 and a 3 are also face-diagonally 
opposite one another in pairs, so a , ai, a 2 , a 3 are the vertices of a regular tetrahedron, T a , of edge length d 
say. In fact, let us choose origin at a , scaling d — \/2 and orient the coordinate axes so that ai = (0, 1, 1), 
&2 = (1,0, 1) and a 3 = (1, 1,0). If we let be the other vertices of the hexahedron, with diagonally 
opposite a i; then the b, also form a tetrahedron, Tt, with edge length d, but of the opposite orientation. So 
in fact we may write hi = b — L&i and L will then be a proper orthogonal matrix. The conditions on the 
bo and L now are that all the faces must be planar, and must not intersect each other except at the edges 
where they meet. Each face has vertices en , hj , afc , hi where {i,j,k,l} are some permutation of {0, 1,2,3}. 
The condition for face a , b 3 , ai, b 2 to be planar is 

A ± : = [ai,b - ia 2 ,b - La 3 ] = [ai,ia 2 ,La 3 ] - [b ,ai,L(a 2 - a 3 )]. (1) 

Equally, the condition for the face b ,a 2 ,bi,a 3 to be planar is 

B x : = [b - bi,b - a 2 ,b - a 3 ] = [Lai,a 2 ,a 3 ] - [b ,iai,a 2 - a 3 ]. (2) 

So the 6 planarity conditions are Ai, Bi and the equations A 2 , A 3 , B 2 , B 3 obtained from them by cyclic 
permutation of the indices {1,2,3}. Since we now have 6 inhomogencous linear equations for the 3 com- 
ponents of bo, consistency of these equations constrains L. It seems computationally easiest to find what 
these constraints are by representing L in terms of a unit quaternion q = qo + qi'i + <j 2 j + g 3 k, so that 
Lx = qxq. If L is rotation by 9 about a unit vector n then q = ±(cos(|#) + nsin(|#)), and we shall choose 
the representation with q > 0. Then a determinant calculation by Mathcmatica shows that consistency of 
equations A\, A 2 , Bi, B 2 gives either qi = ±q 2 or some qi = for i ^ 0. Similarly, by cyclic permutation 
we deduce that either (a) q\ = ±q 2 = ±g 3 , or (b) q\ = or q 2 = or q 3 = 0. Stating this geometrically, in 
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case (a) the axis of the rotation L is parallel to one of the axes of 3-fold rotational symmetry of T a ; while in 
case (b) it is coplanar with one of the pairs of opposite edges of T a . 

In case (a) , if we choose q\ = +q 2 = +93 then the solution is parametrised by q\ , bo = (b, b, b) with 
b = (1 — 4qf)/(l — &q\), b x = b — Lai = b — (4<z 2 , {Qo — <7i) 2 , (90 + 9i) 2 ) an d b 2 , b 3 are obtained from b : by 
cyclic permutation of coordinates. This naturally is the 1-parameter family of solutions constructed earlier 
with 3-fold rotational symmetry about the (1,1,1) direction. The other choices of ± signs give congruent 
hexahedrons to this family. The solution remains valid for — 1/VT2 < q\ < +1/ a/12: at the ends of this 
range when qi = ±l/y/l2 the vertex coincides with hi T i where those suffices are taken from {1,2,3} 
(mod 3), and for \qi\ > l/\/l2 the polyhedron becomes improper because some pairs of faces intersect. An 
example (with qi — 0.2) is shown in Figure 1(a). 




(a) (b) (c) 

In case (a) the dashed line is the axis of 3-fold rotational symmetry. In case (b) the dashed lines mark the 
plane of symmetry. In case (c) the dashed lines mark the 2 planes of symmetry. 

Figure 1: Hexahedrons with all face-diagonals equal. 
In case (b), suppose to be definite we take q% = 0, so q$ + q\ + q\ = 1. Then the coordinates of b are 

bo = (1 + <?o<?2 + <7i<72 - q\ + Zqiqi/lo, 9o - <?o<7i + <7i<?2 + <ll - 1qlq 2 /qo, ql + <7o<7i - 9o<?2 + 2q x q 2 ) . (3) 

This then forms the 2-parameter family constructed earlier, and it can be checked that the face ao, h 2 , a.3, bi 
is congruent to bo,ai,b3,a2, and that there is a plane of symmetry P with b 3 _i being the reflection of a.; 
in P. An example (with q\ — 0.2, q 2 — 0.25) is given in Figure 1(b), where the congruent faces are roughly 
kite-shaped. It remains a valid solution over the region defined by the inequalities 

Wl\ + |®| <j=, (1 - (kil + Id) 2 ) (1 - 2(\ qi \ + \q 2 \r) > 2\ qi q 2 \(\ qi \ + \q 2 \f. (4) 

On the edge of this region P passes through a vertex and so two of the faces fail to be proper quadrilaterals, 
and outside this region the polyhedron becomes improper because some pairs of faces intersect. 

In case (b) when two of the qt vanish, say q 2 — q-3 — 0, then there are two planes of symmetry, and 
the faces are 2 congruent rectangles and 4 congruent symmetric trapezoids, illustrated (for q\ = 0.2) in 
Figure 1(c). 
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B Algorithm for finding a sufficient set of measurements 



In this appendix we give Python code that implements an algorithm to find a sufficient set of measurements. 
The code should be compatible with versions 2.4 and later of Python. The code also uses the 'numpy' 
package, freely available from http://numpy.scipy.org. 

from numpy import array, reshape, concatenate, transpose, compress, \ 

float64, zeros, ones, eye, cross, inner, outer, sometrue 
from numpy. linalg import det, svd 
from math import sqrt, acos, pi, cos, sin 
from itertools import chain, groupby 
from operator import itemgetter 

dcf first (iterable): 

return iterable. next() 

dcf all(conditions): 

for c in conditions: 
if not c: 

return c 
return True 

def any(conditions) : 
for c in conditions: 
if c: 

return c 
return False 

# 

# Numerical linear algebra. We use singular value decomposition to 

# determine the nullspace of a matrix. 

def norm(v): 

return sqrt(inner(v, v)) 

def rank(A, eps = le— 12): 

return sum(svd(A)[l] > le— 12) 

def independent(A, eps = le— 12): 

"""Determine whether the rows of the matrix A are linearly independent""" 
return len(A) == or sum(svd(A)[l] > le-12) == len(A) 

def nullspace(A, eps = le— 12): 

"""null space of the given matrix A 

Find a return a basis for the space of vectors v such that Av = 0.""" 
_, D, VT = svd(A) 

# the second return value (D) of svd is the list of the entries of 

# the diagonal matrix, in decreasing order; if the original matrix 

# has shape (m, n) then D has length min(m, n). We need to pad it 

# out to length n. 

scale = max(D[0], le— 3) 

condition_number = min(a for a in D/scale if a > eps) 



18 



if condition_number <= le— 6: 

# warn that results may be unreliable 
print \ 

"Warning: matrix appears ill— conditioned; " + \ 
"results may be unreliable." + \ 
"Condition number: %s" % condition_number 
print "singular value array: ",D 

return compress(concatenate( 

(D/scale <= eps, [True]*(A.shape[l] - len(D)))), VT, axis=0) 

def independent_subset(labelled_rows, A, eps = le— 12): 

"""find an independent subset from a labelled list of vectors 

Given a labelled list of pairs (label, row) where each row is a 
vector of length n, and an initial m— by— n matrix A of rank m, 
return a list of labels corresponding to a maximal subset of the 
given vectors that, together with the rows of A, forms a linearly 
independent set.""" 

subset = [] 

B = A.copy() 

Bnull = nullspace(B) 

for label, row in labelled_rows: 

if sometrue(abs(inner(Bnull, row)) > eps): 

subset .append(label) 

B = concatenate ((B, [row])) 

Bnull = nullspace(B) 
return subset 

# 

# Planes, lines and points 

# A point in Euclidean 3— space is represented as usual simply as a 

# vector of 3 real numbers; it's also occasionally convenient to think 

# of it as an element of projective 3— space, identifying [x, y, z] 

# with [x : y : z : lj. 

# A line in P3 will be represented simply as a pair of points. 

# The plane with equation ax + by + cz + d = is represented as a 

# quadruple [a, b, c, d] of real numbers; this representation is 

# clearly unique only up to scaling by a nonzero constant. A point 

# [x, y, z] lies on a plane [a, b, c, d] if and only if the inner 

# product of [x, y, z, 1] with [a, b, c, d] vanishes. 

class NoUniqueSolution(Exception) : 

pass 

def plane_containing_points(points): 

"""Find the unique plane containing a given set of points 

Given a list of points in projective 3— space, determine whether they lie 
on a unique plane, and if so return that plane. If not, raise a 
NoUniqueSolution exception." " " 
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v = nullspace(array(points)) 
if len(v) == 1: 

return v[0] 
elif len(v) > 1: 

raise NoUniqueSolution( 

"more than one plane contains the given points") 
elif len(v) < 1: 

raise NoUniqueSolution(" points are not coplanar") 

def line_containing_points (points) : 
v = nullspace(array(points)) 
if len(v) == 2: 

return tuple(v) 
elif len(v) > 2: 

raise NoUniqueSolutionf'more than one line contains the given points") 
elif len(v) < 2: 

raise NoUniqueSolution(" points are not collinear") 

def point _avoiding_planes (planes) : 

"""find a point in P3 not lying on any plane in the given list""" 

n = len(planes) + l 
base = [] 

for plane_set in concatenate( 

(transpose([eye(3, dtype=float64)]*n, [1, 0, 2]), 
— transpose([[range(n)]*3], [1, 2, 0])), 2): 
base.append(first(p for p in plane_set 

if all(independent(base + [p, plane]) for plane in planes))) 
return nullspace(array(base))[0] 



# class for representing permutations of finite sets; this is used for 

# the two permutations of darts required to specify the combinatorial 

# structure of a polyhedron. 

class Permutation(object): 

"""Permutation of a finite set""" 

def __init__(self, pairlist): 

"""initialize from a list of key, value pairs""" 
self.perm_map = dict(pairlist) 
self.domain = set(self.perm_map.keys()) 
if self.domain != set(self.perm_map.values()): 

raise ValueError("map is not a permutation") 

def __eq__ (self, other): 

return self.perm_map == other. perm_map 

def __ne__ (self, other) : 

return not (self == other) 

def __iter__(self): 

return self.perm_map.iteritems() 

def __getitem__(self, key): 

return self.perm_map[key] 
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@staticmethod 
def identity (S): 

"""identity permutation on the given set""" 

return Permutation((s, s) for s in S) 

dcf __muL_(self, other): 

if self.domain == other. domain: 

return Permutation((k, self [ot her [k]]) 

for k in other. perm_map.keys()) 

else: 

raise ValueError(" different domains in multiplication") 

def __str__(self): 

return "" .join(map(str, self.cycles()[0])) 

__repr__ = __str__ 

@staticmethod 

def fromcycle (cycles): 

"""initialize from a list of cycles""" 

return Permutation((c[i— l],c[i]) for c in cycles for i in range (len(c))) 
def cycles (self): 

"""return a list of the cycles of the given permutation, together with 
a function that maps every element of the domain to the cycle 
containing it.""" 

p = diet (self. perm_map) 
def cycle_starting_at(p, v): 
cycle = [] 
while v in p: 

cycle, append (v) 
v = p.pop(v) 
return tuple (cycle) 
cycleJist = [] 
while p: 

cycle _list . append (cycle_st arting_at (p , min (p) ) ) 

cycle_containing = dict((e, c) for c in cycleJist for e in c) 
return cycleJist, (lambda x: cycle_containing[x]) 

def equivalence_classes(S, pairs): 

"""given a set or list S, and a sequence of pairs (s, t) of elements of S, 
compute the equivalence classes of S generated by the relations s ~ t for 
(s, t) in pairs. 

Returns a list of equivalence classes, with each equivalence class 
given as a list of elements of S. 

Example: 

>>> equivalence_classes(range(10), [(1, 4), (2, 6), (2, 3), 

(0, 7), (1, 8), (3, 0)]) 

[[0, 2, 3, 6, 7], [1, 4, 8], [5], [9]] 
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?! !) 5? 

# Algorithm based on that described in 'Numerical Recipes in C, 

# 2nd edition', section 8.6; in turn based on Knuth. 

def root(s): 

while s != parent [s]: 

s = parent [s] 
return s 

parent = dict((s, s) for s in S) 
for s, t in pairs: 

s, t = root(s), root(t) 

if s != t: 

parent[max(s, t)] = min(s, t) 

# by this stage, two elements s and t are in the same equivalence 

# class if and only if they have the same root. So we just 

# have to collect together items with the same root. 

return [map(itemgetter(l), c[l]) for c in 

groupby(sorted((root(s), s) for s in S), itemgetter(O))] 

# 

class CombinatorialError(Exception) : 

def __init__(self, message, abstract_poly) : 
self, message = message 
self, abstract _poly = abstract_poly 

def __str__(self, message, abstract_poly): 
return message 

class AbstractPolyhedron(object) : 

"""Combinatorial data associated to an (oriented) polyhedron. 

The fundamental element of the combinatorial information is the 
'dart'. One dart is associated to each pair (v, f) where v is a 
vertex, f is a face, and v is incident with f. Thus each dart has 
an associated vertex and face, each vertex can be described by the 
set of darts around it, and each face can be described by the set 
of darts on it. Each dart is also naturally associated to an 
edge: in general for a dart d = (v, f) there will be two edges 
incident with both v and f; let w be the next vertex from v on f, 
moving *counterclockwise* around f, then we associate the edge v— w 
to the dart (v, f). 

Now the combinatorial data needed are completely described by two 
permutations: one that associates to each dart d the next dart 
counterclockwise that's associated to the same vertex, and one 
that associates to each dart d the next dart counterclockwise 
that's associated to the same face. 

Attributes of an AbstractPolyhedron instance: 

darts : a list of the underlying darts of the polyhedron; these can be 
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instances of any class that's equipped with an ordering, for example 

integers or strings, 
vertexcycle : a Permutation instance such that for each dart d, 

associated to a vertex v, vertexcycle [d] is the next dart 

around v from d, moving counterclockwise around v. 
facecycle: similar to vertexcycle: gives the permutation of darts 

around each face, again counterclockwise, 
edgecycle: composition of vertexcycle with facecycle. Sends 

a dart d to the other dart associated to the same edge. 

vertices: list of vertices of the polyhedron. Each vertex is 

(currently) represented as a list of darts, 
faces: list of faces of the polyhedron 
edges: list of edges of the polyhedron 

components : list of components of the polyhedron; each component 
is represented as a list of darts. For a connected (nonempty) 
polyhedron, there will be only one component. 

euler : the Euler characteristic of the abstract polyhedron 

V : function that associates to a given dart d the corresponding 
vertex 

F : function associating to each dart d the corresponding face 
E : function associating to each dart d the corresponding edge 

55 55 55 

def __init__(self, vertexcycle, facecycle): 

"""given permutations 'vertexcycle' and 'facecycle' describing 
the vertex and face permutations on a set of darts, validate 
and construct the corresponding abstract (oriented) polyhedron. 

55 55 55 

if vertexcycle. domain == facecycle. domain: 
self.darts = vertexcycle. domain 

else: 

raise CombinatorialError("the vertex and face permutations "\ 
"should have the same underlying set of darts", self) 

self, edgecycle = vertexcycle * facecycle 
self, facecycle = facecycle 
self, vertexcycle = vertexcycle 

if self.edgecycle*self.edgecycle != Permutation. identity (self.darts): 

raise CombinatorialError(" every cycle of the edge permutation "\ 
"should have length <= 2", self) 

self. vertices, self.V = self. vertexcycle. cycles() 
self.faces, self.F = self.facecycle.cycles() 
self.edges, self.E = self.edgecycle.cycles() 
self, components = equivalence_classes(self.darts, 
chain(self. vertexcycle, self. facecycle)) 

self.euler = 2 — len(self.darts) — 2*len(self.components) \ 
+ len(self. vertices) + len(self.edges) + len(self.faces) 
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# additional checks that apply to a connected polyhedron, 

# embeddable in R3 with planar faces. 
self.validate() 

# some basic combinatorial operations; others can be defined similarly 
def vertices_on_face(self, f): return map(self.V, f) 

def vertices_on_edge(self, e): return map(self.V, e) 
def faces_around_vertex(self, v): return map(self.F, v) 
def faces_around_edge(self, e): return map(self.F, e) 
def edges_on_face(self, f): return map(self.E, f) 
def edges_around_vertex(self, v): return map(self.E, v) 

def __str__(self): 
return """\ 
Abstract polyhedron with: 
dart set : %s 
vertex permutation : %s 
face permutation : %s 
edge permutation : %s 
""" % (self.darts, self.vertexcycle, self.facecycle, self.edgecycle) 

__repr__ = __str__ 

def validate (self): 

# each edge should contain at least two darts; each vertex and face 

# should have at least three. 

bad_edge = any(len(c) < 2 for c in self.edges) 
if bad_edge: 

raise CombinatorialError("edge %s has the wrong length: it "\ 
"should involve exactly two darts" % bad_edge, self) 

bad_vertex = any(len(c) < 3 for c in self. vertices) 
if bad_vertex: 

raise CombinatorialError(" vertex %s has the wrong length: it "\ 
"should involve at least three darts" % bad_vertex, self) 

bad _face = any(len(c) < 3 for c in self.faces) 
if bad_face: 

raise CombinatorialError("face %s has the wrong length: it "\ 
"should involve at least three darts" % bad_face, self) 

# verify that the abstract polyhedron is connected 
if len(self.components) != 1: 

raise CombinatorialError(" Disconnected polyhedron. " \ 
" Components: %s" % self.components, self) 

class Polyhedron(AbstractPolyhedron) : 

# a polyhedron should consist of: 

# (1) an underlying abstract polyhedron 

# (2) an embedding of that abstract polyhedron 

def __init__(self, vertexcycle, facecycle, embedding, allowScaling = False): 

# initialize the combinatorial data 

AbstractPolyhedron.__init__(self, vertexcycle, facecycle) 
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# store position information for vertices (in PS) 
vertexeqn = dict( 

(self.V(d), array(pos + [1.])) 
for d, pos in embedding. items()) 

# verify that every vertex has had its position given 

# exactly once 

if not (len(vertexeqn) == len(embedding) == len(self. vertices)): 
raise ValueError(" missing vertices or repeated vertex "\ 
"in given embedding") 

# find the coordinates of the plane for each face (in P3) 

# this raises an error if the set of given points is not coplanar 
faceeqn = {} 

for f in self. faces: 

vertex_positions = [vertexeqn[v] for v in self.vertices_on_face(f)] 
try: 

faceeqn[f] = plane_containing_points(vertex_positions) 
except NoUniqueSolution: 

raise ValueError("the points of face %s do not lie " 
"in a unique plane" % (str(f))) 

# choose a reference point in P3 that doesn't lie on any of the planes 

# and doesn't lie on the plane at infinity 

basept = point_avoiding_planes(faceeqn.values() + [array([0., 0., 0., 1.])]) 

# project into R3, shifting so that the reference point goes 

# to (0, 0, 0) Then none of the faces goes through 0, 0, 0, so 

# each face can be represented by a triple a, b, c, 

# corresponding to the equation ax + by + cz + 1 = 0. 
vertexpos = dict((v, p[:3]/p[3] — basept[:3]/basept[3]) 

for v, p in vertexeqn. itemsQ) 
facepos = dict((f, e[:3]*basept[3]/inner(e, basept)) 
for f, e in faceeqn. items()) 

# now fill in the relationmatrix (= D(phi)) 
ambient_dim = 3*(len(self.vertices)+len(self.faces)) 
relationmatrix = zeros ((len (self. darts), ambient_dim), dtype = float64) 
for i, d in enumerate(self.darts): 

# construct the row corresponding to a vertex face pair 

v = self.V(d) 
f = self.F(d) 

vi = 3*self. vertices. index(v) 

fi = 3*(len(self.vertices)+self.faces.index(f)) 
relationmatrix[i, vi:vi+3] = facepos[f] 
relationmatrix[i, fi:fi+3] = vertexpos[v] 

# check to see whether the relation matrix has full rank 
if rank (relationmatrix) != len(relationmatrix): 

print ("Warning: relation matrix has rank %d; " 
"expected rank %d" % 

(rank(relationmatrix), len(relationmatrix))) 

# compute the vectors of the Lie group: shifts, rotations and scalings 
shifts = [eye(3, dtype=float64) for v in self. vertices] + \ 
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[outer(facepos[f], facepos[f|) for f in self.faces] 
shifts = reshape(transpose(shifts, [1, 0, 2]), (3, ambient_dim)) 
rotations = \ 

[cross(eye(3, dtype=fioat64) , vertexpos[v]) for v in self, vert ices] + \ 
[cross(eye(3, dtype=fioat64) , facepos[f]) for fin self.faces] 

rotations = reshape(transpose(rotations, [1, 0, 2]), (3, ambient _dim)) 

if allowScaling: 

scalings = [vertexpos[v] for v in self. vertices] + \ 

[— facepos[f] for f in self.faces] 
scalings = reshape (scalings, (1, ambient _dim)) 
liealgebra = concatenate((shifts, rotations, scalings)) 

else: 

liealgebra = concatenate((shifts, rotations)) 

self.euler_characteristic = \ 

len(self.faces) + len(self. vertices) — len(self.edges) 

# Polyhedron attributes: 

# darts : list of darts ( = vertex— face pairs ) 

# vertices : list of darts representing the vertices, one per vertex 

# faces : list of darts representing the faces, one per face 

# liealgebra contains the vectors of the lie algebra 

# vertexcyclefd] : the dart obtained by moving counterclockwise 

# about the vertex associated to d 

# facecyclefd] : the dart obtained by moving counterclockwise 

# about the face associated to d 

# vertexposfv] : the position of vertex v in R3 ( after shifting to 

# make sure that no faces contain 0, 0, 0) 

# faceposff] : the plane containing face f (after shifting), given 

# by a triple (a, b, c) representing the plane with equation 

# ax + by + cz + 1 = 0. 

# embeddingfv] : the original position of vertex v 

# faceeqnff] : the equation of face f (before shifting), given 

# by a 4~ tuple (a, b, c, d) representing the plane with 

# equation ax + by + cz + d = 0. 

self, liealgebra = liealgebra 
self.vertexpos = vertexpos 
self.facepos = facepos 
self. allowScaling = allowScaling 
self.ambient_dim = ambient_dim 
self.relationmatrix = relationmatrix 
self.vertexeqn = vertexeqn 
self.faceeqn = faceeqn 

dcf filter _measurements(self, measurementlist): 

measurements_used = independent_subset([(m, m.derivative()) 
for m in measurementlist], self.relationmatrix) 

defect = self, ambient _dim — len(self.liealgebra) — \ 

len(self.relationmatrix) — len(measurements_used) 
return measurements_used, defect 

def length(self, vl, v2): 

return LengthMeasurement(self, vl, v2) 
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def angle(self, vl, v2, v3): 

return AngleMeasurement(self, vl, v2, v3) 

def dihedral(self, fl, f2): 

return DihedralMeasurement(self, fl, f2) 

class Measurement(object): 

"""Class representing a particular measurement of a 
polyhedron.""" 

# we '11 subclass this for length, angle and dihedral angle 

# measurements 

# each subclass should implement value(), derivativef), str() 
pass 

# Note that the length measurement actually only needs to depend on 

# the abstract polyhedron; we need a genuine embedded polyhedron only 

# to apply the measurement. 

class LengthMeasurement(Measurement) : 

"""Class representing the measurement of the distance 
between two vertices of a polyhedron.""" 

def __init__(self, P, vl, v2): 

if not (vl in P. vertices and v2 in P.vertices): 
raise ValueError(" invalid vertices") 

if P.allowScaling: 

raise ValueError("working up to similarity only; " 

"length measurements not permitted") 

# P is the polyhedron; vl and v2 are the relevant vertices 
self.P = P 

self.vl, self.v2 = vl, v2 
dcf __str__(self): 

return "Length measurement V(%s) V(%s)" % (self.vl [0], self.v2[0]) 

def value (self): 

return norm(self.P.vertexpos[self.vl] — self.P. vertexpos[self.v2]) 

def derivative (self): 
P = self.P 

pi, p2 = P.vertexpos[self.vl], P.vertexpos[self.v2] 
cliff = pi - p2 
length = norm(diff) 

row = zeros((len(P. vertices) + len(P. faces), 3), dtype=float64) 
row[P. vertices. index(self.vl)] += diff/length 
row[P. vertices. index(self.v2)] — = diff/length 
return reshape(row, (P.ambient_dfm,)) 

class AngleMeasurement (Measurement): 



27 



# cosine of the angle v2—vl—v3 
dcf __init__(self, P, vl, v2, v3): 

if not all(v in P. vertices for v in [vl, v2, v3]): 

raise ValueError(" invalid vertices") 
self.P = P 

self.vl, self.v2, self.v3 = vl, v2, v3 
def __str__(self): 

return "Angle measurement V(%s) V(%s) V(%s)" % \ 

(self.v2[0], self.vl [0], self.v3[0]) 

def value (self): 
P = self.P 

(pi, p2, p3) = (P.vertexpos[self.vl], 

P.vertexpos[self.v2], P.vertexpos[self.v3]) 
diff2, diff3 = p2 - pi, p3 - pi 
return inner(diff2, diff3)/norm(diff2)/norm(diff3) 

def derivative (self): 

"""vector representing the derivative of the cosine of the angle 
measurement between vectors v2— vl and v3— vl""" 

P = self.P 

(pi, p2, p3) = (P.vertexpos[self.vl], 

P.vertexpos[self.v2], P.vertexpos[self.v3]) 
diff2, diff3 = p2 - pi, p3 - pi 
a23 = inner(diff2, diff3) 
a22 = inner(diff2, diff2) 
a33 = inner(diff3, diff3) 
M = sqrt(a22 * a33) 

cosangle = a23/M 

dvl = (diff2*a23/a22 + diff3*a23/a33 - diff2 - diff3)/M 
dv2 = (diff3 - diff2*a23/a22)/M 
dv3 = (diff2 - diff3*a23/a33)/M 

row = zeros((len(P. vertices) + len(P. faces), 3), dtype = float64) 
row[P. vertices. index(self.vl)] += dvl 
row[P. vertices. index(self.v2)] += dv2 
row[P. vertices. index(self.v3)] += dv3 
return reshape(row, (P.ambient_dim,)) 

class DihedralMeasurement(Measurement) : 

# measurement of the angle between faces fl and f2 

dcf __init__(self, P, fl, f2): 

if not all(f in P.faces for f in [fl, f2]): 
raise ValueError(" invalid faces") 
self.P = P 

self.fl, self.f2 = fl, f2 
def __str__(self): 

return "Dihedral measurement F(%s) F(%s)" % (self.fl [0], self.f2[0]) 

def value (self): 
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# angle between two faces: the faces have corresponding embeddings 
nl, n2 = self.P.facepos[self.fl], self.P.facepos[self.f2] 



# nl and n2 are (not necessarily unit) normal vectors to the 

# two faces note that nl and n2 don 't necessarily have 

# corresponding orientations: one of nl and n2 might point 

# outwards while the other points inwards. We need to fix 

# this. For now we just compute the (cosine of the) angle 

# between nl and n2. 

return inner(nl, n2)/norm(nl)/norm(n2) 

dcf derivative (self): 
P = self.P 

nl, n2 = P.facepos[self.fl], P.facepos[self.f2] 

all = inner(nl, nl) 
al2 = inner(nl, n2) 
a22 = inner(n2, n2) 

dnl = (n2 - nl*al2/all)/sqrt(all*a22) 
dn2 = (nl - n2*al2/a22)/sqrt(all*a22) 

row = zeros((len(P. vertices) + len(P. faces), 3), dtype=fioat64) 
row[len(P. vertices) + P.faces.index(self.fl)] += dnl 
row[len(P. vertices) + P.faces.index(self.f2)] += dn2 
return reshape(row, (P.ambient_dim,)) 

# 

# Output routines 

def print_polyhedron(P): 

print "Vertices (number, position, cycle)\n" \ 



print "".join("%2d " % v[0] + 
" (% 7.4f, % 7.4f, % 7.4f) " % 

tuple(P.vertexeqn[v][:3]/P.vertexeqn[v] [3]) + 
str(v) + "\n" for v in P. vertices) 

print "Faces (number, equation, cycle)\n" \ 



print "".join("%2d " % f[0] + 

"% 7.4fx + % 7.4fy + % 7 Ah = % 7 At " % tuple(P.faceeqn[f]) + 
str(f) + "\n" for f in P. faces) 

# print results in a meaningful format 
def print_results(P, results): 

for r in results: 

print "%s: expected value: % 7.4P \ 
% (r, r.valueQ) 

# 

# Platonic solids 

tetrahedron = Polyhedron( 

vertexcycle = Permutation. fromcycle([ 
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[0, 6 ,3], [1, 5, 9], [2, 11, 7], [4, 8, 10]]), 
facecycle = Permutation. fromcycle( [range (i, i+3) for i in range(0, 12, 3)]), 
embedding = dict([ 

[0, [1., 1, 1.]], 

[1, [1., -1., -1.]], 

[2, [-1, L, -1.]], 

[4, hi-, -1., 1.]]])) 

cube = Polyhedron( 

vertexcycle = Permutation. fromcycle([ 

[0, 4, 8], [1, 11, 22], [2, 21, 19], [3, 18, 5], 

[6, 17, 15], [7, 14, 9], [10, 13, 23], [12, 16, 20]]), 
facecycle = Permutation. fromcycle( 

[range(i, i+4) for i in range(0, 24, 4)]), 
embedding = dict([ 

[0, [1., 1, 1.]], 

[1, [1, 1., -1.]], 

[2, [1., -1., -1.]], 

[3, [1, -1, 1.]], 

[6, [-1, -1., 1.]], 

[7, [-1., 1, 1.]], 

[10, [-1., 1, -1.]], 

[12, [-I., -1., -1.]]])) 

twisted_cube = Polyhedron( 

vertexcycle = Permutation. fromcycleQ 

[0, 4, 8], [1, 11, 22], [2, 21, 19], [3, 18, 5], 

[6, 17, 15], [7, 14, 9], [10, 13, 23], [12, 16, 20]]), 
facecycle = Permutation. fromcycle( 

[range(i, i+4) for i in range(0, 24, 4)]), 
embedding = dict([ 

[0, [59./51., 451./510., 451./510.]], 

[1, [1., 1, 0.]], 

[2, [13./17., -13./170., -13./170.]], 

[3, [1., 0., 1.]], 

[6, [0., .1, 1.1]], 

[7, [0., 1, 1.]], 

[10, [0., 1.1, 0.1]], 

[12, [0., 0., 0.]]])) 

twisted_cube2 = Polyhedron( 

vertexcycle = Permutation. fromcycleQ 

[0, 4, 8], [1, 11, 22], [2, 21, 19], [3, 18, 5], 

[6, 17, 15], [7, 14, 9], [10, 13, 23], [12, 16, 20]]), 
facecycle = Permutation. fromcycle( 

[rangc(i, i+4) for i in range(0, 24, 4)]), 
embedding = dict([ 

[0, [-1.1, -1.1, 1.]], 

[1, [1, 1, -1.]], 

[2, [-1.1, .9, -1.]], 

[3, [1., -1, 1.]], 

[6, [.9, .9, 1.]], 

[7, [-1, L, 1.]], 

[10, [.9, -1.1, -1.]], 

[12, [-1., -1., -1.]]])) 
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# rotate 2nd tetrahedron by 45 degrees about z—axis 

# transformation: [x, y, z] — > [x cost — y smt, x sint + y cost, zj 
theta = 0.3 

cost = cos(theta) 
sint = sin(theta) 
twisted_cube3 = Polyhedron( 

vertexcycle = Permutation. fromcycleQ 

[0, 4, 8], [1, 11, 22], [2, 21, 19], [3, 18, 5], 

[6, 17, 15], [7, 14, 9], [10, 13, 23], [12, 16, 20]]), 
facecycle = Permutation. fromcycle( 

[rangc(i, i+4) for i in range(0, 24, 4)]), 
embedding = dict([ 

[0, [1, 1., 1.]], 

[2, [1, -1., -1.]], 

[6, [-1., -1., 1.]], 

[10, [-1., 1, -1.]], 

[1, [cost— sint, sint+cost, — 1.]], 

[3, [cost+sint, sint— cost, 1.]], 

[7, [—cost— sint, —sint+cost, 1.]], 

[12, [—cost+sint, — sint— cost,— 1.]]])) 

octahedron = Polyhedron( 

vertexcycle = Permutation. fromcycle( 

[range(i, i+4) for i in range(0, 24, 4)]), 
facecycle = Permutation. fromcycle([ 

[0, 4, 8], [1, 11, 22], [2, 21, 19], [3, 18, 5], 

[6, 17, 15], [7, 14, 9], [10, 13, 23], [12, 16, 20]]), 
embedding = dict([ 

[0, [ 1, 0., 0.]], 

[4, [0., 0., 1.]], 

[8, [0., L, 0.]], 

[12, [-1., 0., 0.]], 

[16, [0., -1., 0.]], 

[20, [0., 0., -1.]]])) 

t, tr = (sqrt(5) + l.)/2., (sqrt(5)-l.)/2. 
dodecahedron = Polyhedron( 

vertexcycle = Permutation. fromcycleQ 

[0, 5, 11], [1, 10, 16], [2, 15, 21], [3, 20, 26], [4, 25, 6], 

[7, 29, 58], [8, 57, 54], [9, 53, 12], [13, 52, 49], [14, 48, 17], 

[18, 47, 44], [19, 43, 22], [23, 42, 39], [24, 38, 27], [28, 37, 59], 

[30, 35, 41], [31, 40, 46], [32, 45, 51], [33, 50, 56], [34, 55, 36]]), 
facecycle = Permutation. fromcycle( 

[rangc(i, i+5) for i in range(0, 60, 5)]), 
embedding = dict([ 

[0, [1., 1., 1.]], 

[1, [t, tr, 0.]], 

[2, [1., 1., -1.]], 

[3, [0., t, -tr]], 

[4, [0., t, tr]], 

[7, [-1., 1., 1.]], 

[8, [-tr, 0., t]], 

[9, [tr, 0., t]], 

[13, [1., -1., 1.]], 

[14, [t, -tr, 0.]], 

[18, [1., -1., -1.]], 
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[19, [tr, 0., -t]], 
[23, [-tr, 0., -t]], 
[24, [-1., 1, -1.]], 
[28, [-t, tr, 0.]], 
[30, [-1, -1, -1.]], 
[31, [0., -t, -tr]], 
[32, [0., -t, tr]], 
[33, [-1., -1, 1.]], 
[34, [-t, -tr, 0.]]]), 
allowScaling = True) 

icosahedron = Polyhedron( 

vertexcycle = Permutation. fromcycle( 

[range(i, i+5) for i in range(0, 60, 5)]), 
facecycle = Permutation. fromcycle([ 

[0, 5, 11], [1, 10, 16], [2, 15', 21], [3, 20, 26], [4, 25, 6], 

[7, 29, 58], [8, 57, 54], [9, 53, 12], [13, 52, 49], [14, 48, 17], 

[18, 47, 44], [19, 43, 22], [23, 42, 39], [24, 38, 27], [28, 37, 59], 

[30, 35, 41], [31, 40, 46], [32, 45, 51], [33, 50, 56], [34, 55, 36]]), 
embedding = dict([ 

[0, [1., t, 0.]], 

[5, [0., 1, t]], 

[10, [t, 0., 1.]], 

[15, [t, 0., -1.]], 

[20, [0., 1, -t]], 

[25, [-1, t, 0.]], 

[30, [-1, -t, 0.]], 

[35, [-t, 0., -1.]], 

[40, [0., -1, -t]], 

[45, [1, -t, 0.]], 

[50, [0., -1., t]], 

[55, [-t, 0., 1.]]])) 



# a tetrahedron with a notch taken out of the center of one edge. The 

# result is a polyhedron with 12 edges, 6 faces and 8 vertices; the 

# abstract polyhedron can be realized as a simple polyhedron, but not 

# as a convex polyhedron ( the corresponding planar graph is not 

# 3— connected, so Steinitz' theorem doesn't apply). 

notched_tetrahedron = Polyhedron( 

vertexcycle = Permutation. fromcycleQ 

[0, 12, 6], [1, 11, 18], [2, 20, 22], [3, 21, 9], 

[4, 8, 17], [5, 16, 13], [7, 14, 15], [10, 23, 19]]), 
facecycle = Permutation. fromcycle([ 

[0, 1, 2, 3, 4, 5], [6, 7, 8, 9, 10, 11], 

[12, 13, 14], [15, 16, 17], [18, 19, 20], [21, 22, 23]]), 
embedding = dict([ 

[0, [1., 1., 1.]], 

[1, [0.5, 1., 0.5]], 

[2, [-0.5, 0., 0.5]], 

[3, [-0.5, 1., -0.5]], 

[4, [-1., 1., -1.]], 

[5, [-1., -1., 1.]], 

[7, [1., -1., -1.]], 

[10, [0.5, 0., -0.5]]])) 
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square_pyramid = Polyhedron( 

vertexcycle = Permutation. fromcycle([ 

[0, 3, 6, 9], [1, 11, 13], [2, 12, 4], [5, 15, 7], [8, 14, 10]]), 
facecycle = Permutation. fromcycle([ 

[0, 1, 2], [3, 4, 5], [6, 7, 8], [9, 10, 11], [12, 13, 14, 15]]), 
embedding = dict([ 

[0, [0., 0., 1.]], 

[1, [1., 1., 0.]], 

[2, [-1., 1., 0.]], 

[5, [-1., -1., 0.]], 

[8, [1., -1., 0.]]])) 

# square pyramid with two of the triangular faces equal 

degenerate_square_pyramid = Polyhedron( 
vertexcycle = Permutation. fromcycle([ 

[0, 3, 6, 9], [1, 11, 13], [2, 12, 4], [5, 15, 7], [8, 14, 10]]), 
facecycle = Permutation. fromcycle([ 

[0, 1, 2], [3, 4, 5], [6, 7, 8], [9, 10, 11], [12, 13, 14, 15]]), 
embedding = dict([ 

[0, [0., 0., 1.]], 

[1, [1., 1., 0.]], 

[2, [0., 0., 0.]], 

[5, [-1., -1., 0.]], 

[8, [1., -1., 0.]]])) 

# torus (embedding as polyhedron with 9 quadrilateral faces) defined 

# for experimentation purposes only; since it's not a convex 

# polyhedron, the assumptions of the paper and the algorithm don't 

# apply. Expect the PolyhedronQ call for the torus to throw up a 

# warning: "Warning: relation matrix has rank 35; expected rank 36". 

s3 = sqrt(3.) 

hs3 = s3/2. 

torus = Polyhedron( 

facecycle = Permutation. fromcycle( 

[range(i, i+4) for i in range(0, 36, 4)]), 
vertexcycle = Permutation. fromcycle([ 

[0, 9, 34, 27], 

[1, 26, 31, 4], 

[2, 7, 16, 13], 

[3, 12, 21, 10], 

[5, 30, 35, 8], 

[6, 11, 20, 17], 

[14, 19, 28, 25], 

[15, 24, 33, 22], 

[18, 23, 32, 29]]), 
embedding = dict([ 

[0, [1., 0., 0.]], 

[1, [-.5, hs3, 0.]], 

[2, [-1., S3, 1.]], 

[3, [2., 0., 1.]], 

[5, [-.5, -hs3, 0.]], 

[6, [-1., -S3, 1.]], 

[14, [-1., s3, -1.]], 
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[15, [2., 0., -1.]], 

[18, [-1., -s3, -1.]]])) 

# szilassi polyhedron; combinatorial information based on the heawood graph 

szilassi = AbstractPolyhedron( 

vertexcycle = Permutation. fromcycle( 

[range(i, i+3) for i in range(0, 42, 3)]), 
facecycle = Permutation. fromcycle( 

[f(i + j) % 42 for i in [0, 5, 31, 27, 26, 40]] for j in range(0, 42, 6)])) 

# 

# attempt to find a set of angles that's sufficient for a dodecahedron; 

# we can do this using face measurements only 

D = dodecahedron 
faceangles = [] 
for f in D. faces: 

vs = D.vertices_on_face(f) 
for i in range(5): 

for j in range(i+l, 5): 

for k in range (j+1, 5): 

faceangles. append(AngleMeasurement(D, vs[i], vs[j], vs[k])) 
faceangles. append(AngleMeasurement(D, vs[j], vs[k], vs[i])) 

results, defect = dodecahedron. filter_measurements(faceangles) 
print "Regular dodecahedron" 
print _polyhedron(dodecahedron) 
if defect == 0: 

print "A sufficient set of angle measurements for the dodecahedron" 
print_results(dodecahedron, results) 

else: 

print ("The following measurements are independent, " 

"but not infmitesimally sufficient.") 
print "%d more measurements are required" % defect 
print_results(dodecahedron, results) 
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